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Battery  design  is  a  critical  aspect  of  material  and  system  development  that  leads  to  the  commercializa¬ 
tion  of  effective  electrochemical  energy  storage  systems.  Successful  modeling  of  battery  designs  relies 
upon  accurate  calculation  of  the  area  specific  impedance  (ASI).  A  simplified  calculation  of  the  ASI  is 
presented  that  accounts  for  physical  limitations  without  performing  computationally  expensive  calcu¬ 
lations.  The  limiting  currents  for  transport  within  the  electrolyte  and  within  the  intercalation  materials 
are  implemented  into  a  linear  form  of  the  Butler-Volmer  equation  to  calculate  the  interfacial  impedance. 
Lithium-ion  batteries  are  then  designed  to  examine  the  effect  of  power  to  energy  ratio  on  battery  dimen¬ 
sions.  A  large  ASI  is  shown  to  be  detrimental  to  battery  design  regardless  if  the  increase  in  impedance 
results  from  mass  transport  limitations  or  a  reduction  in  electrochemical  active  area  due  to  small  elec¬ 
trode  loadings.  The  smaller  electrochemical  active  area  does  not  increase  the  voltage  losses  of  a  battery 
when  a  constant  C-rate  is  maintained.  However,  the  higher  ASI  values  from  low  electrode  loadings  require 
a  larger  separator  and  current  collector  area  resulting  in  a  greater  battery  volume  and  weight  to  achieve 
similar  energy  and  power  requirements  when  compared  to  a  system  with  a  lower  ASI. 
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1.  Introduction 

The  ability  to  calculate  and  project  battery  weight,  volume, 
material  requirements,  and  cost  is  extremely  valuable  to  the 
research,  industrial,  and  economic  communities  alike.  This  calcu¬ 
lation  is  even  more  valuable  if  it  is  computed  in  an  efficient  and 
time  effective  manner  while  still  producing  meaningful  and  reason¬ 
ably  accurate  results.  To  achieve  this  end,  we  have  been  developing 
relatively  simple  models  that  use  experimentally  measured  elec¬ 
trochemical  data,  physical  properties  and  proposed  cell  designs  as 
inputs  to  the  calculation.  The  model  then  allows  the  user  to  evalu¬ 
ate  the  battery  design  for  various  cell  chemistries  to  meet  various 
power  and  energy  requirements. 

The  required  experimental  data  is  acquired  through  straight¬ 
forward  experiments.  The  reversible  capacity  and  voltage  profile 
of  the  positive  and  negative  electrode  materials  at  the  cycling  rate 
of  interest  (e.g.  C  or  C/3)  must  be  measured  but  is  often  available 
in  the  literature.  Additionally,  measurements  of  the  area  specific 
impedance  (ASI)  for  the  discharge  rate  used  for  energy  and  the  rate 
used  for  power  are  required.  In  transportation  applications,  the  bat¬ 
tery  must  be  able  to  achieve  specified  pulse  discharge  and  charge 
conditions.  The  ASI  for  pulse  power  transportation  applications  can 
be  measured  using  the  hybrid  pulse  power  characterization  (HPPC) 
technique  [1  ].  An  example  of  the  current  and  voltage  behavior  dur¬ 
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ing  a  HPPC  test  may  be  found  in  Fig.  1.  In  the  HPPC  technique,  the 
ASI  is  measured  as  a  function  of  depth  of  discharge  at  10%  inter¬ 
vals.  The  current  is  pulsed  in  discharge  mode  for  a  period  of  time, 
presently  10  s.  Then  the  current  is  stopped  and  voltage  is  allowed 
to  relax.  Next,  the  current  is  pulsed  in  the  charge  mode  at  typi¬ 
cally  75%  of  the  discharge  value.  The  ASI  is  calculated  separately  for 
charge  and  discharge  pulses,  Eq.  (1).  V  stands  for  the  cell  voltage 
and  I  represents  the  current  density.  The  subscripts  refer  to  value 
before  the  pulse,  tO,  and  at  the  end  of  the  pulse,  tl. 


The  ASI  is  not  an  inherent  constant  of  a  specific  battery  chem¬ 
istry  or  cell  design.  The  measured  value  of  the  ASI  is  a  complex 
combination  of  resistances  within  the  battery  resulting  from  the 
physical  processes  occurring  at  different  length  and  time  scales. 
Therefore,  the  measured  value  is  a  function  of  many  factors  (state 
of  charge,  pulse  length,  current  density,  C-rate,  particle  size,  trans¬ 
port  and  kinetic  parameters,  etc).  Three  general  approaches  exist 
to  account  for  the  ASI  in  battery  design  modeling.  First,  one  could 
systematically  design  and  assemble  a  large  number  of  cells  with 
changing  parameters  and  measure  the  resulting  ASIs.  Secondly,  a 
complex  physics  based  model  could  be  used  to  predict  the  ASI  for 
different  electrode  designs  and  then  this  could  be  used  to  predict 
a  “scaled-up”  battery  design  [2,3].  Finally,  one  could  take  the  ASI 
to  be  constant  over  the  range  of  designs  studied.  The  last  option 
is  the  simplest  and  thus  is  most  commonly  used  in  battery  design 
models.  Our  previous  battery  design  models  assumed  the  ASI  for 


0378-7753 /$  -  see  front  matter  ©  2010  Elsevier  B.V.  All  rights  reserved, 
doi:  10.1 01 6/j.jpowsour.201 0.1 0.020 


2290 


ICG.  Gallagher  et  al.  /  Journal  of  Power  Sources  196  (201 1 )  2289-2297 


Nomenclature 

a 

electrochemically  active  area  per  unit  electrode  vol¬ 
ume,  cm-1 

A 

area,  cm2 

ASI 

area  specific  impedance,  £2  cm2 

9 

concentration  of  lithium  in  region  j,  mol  cm-3 

Cj 

reference  concentration  of  lithium  in  region  j, 
mol  cm-3 

c 

cell  capacity,  Ah 

DUP¥6 

diffusion  coefficient  of  LiPF6  based  electrolyte, 
cm2  s-1 

Ds 

diffusion  coefficient  of  lithium  in  active  material 
particle,  cm2  s-1 

F 

Faraday’s  constant,  96485.3  C  mol-1 

io 

exchange  current  per  unit  of  interfacial  area,  A  cm-2 

h 

superficial  solution  current  in  electrolyte  phase, 

A  cm-2 

in 

transfer  current  per  unit  of  interfacial  area,  A  cm-2 

I 

current  density,  A  cm2 

rionic 

him 

limiting  current  density  for  lithium  ion  transport 
through  separator,  A  cm2 

rssDiff 

him 

limiting  current  density  for  lithium  diffusion  in 
active  material  particles,  A  cm2 

L 

region  thickness,  cm 

^cell 

number  of  cells  in  the  battery 

P 

battery  power,  W 

rc 

C-rate,  h-1 

rP 

diffusion  length  of  active  material,  cm 

R 

universal  gas  constant,  8.31 44  J  mol-1  K-1 

fpulse 

duration  of  current  pulse,  s 

t+ 

transference  of  lithium  ions  with  respect  to  the  sol¬ 
vent 

S 

specific  surface  area,  m2  g-1 

T 

absolute  temperature,  K 

a 

specific  capacity,  Ah  g-1 

u 

equilibrium  potential,  V 

V 

voltage,  V 

y 

molar  fraction  of  lithium  in  solid  particles 

z+ 

charge  of  lithium  ion 

aa 

anodic  transfer  coefficient 

ac 

cathodic  transfer  coefficient 

S 

thickness,  cm 

£ 

volume  fraction 

v 

surface  overpotential,  V 

K 

conductivity,  S  cm-1 

P 

active  material  density,  g  cm-3 

Subscript  and  superscript 

act 

active  material 

batt 

battery 

const 

constant 

inact 

inactive  material 

intf 

interfacial 

ocv 

open  circuit  voltage 

pos 

positive  electrode 

neg 

negative  electrode 

sep 

separator 

to 

before  the  pulse 

tl 

at  the  end  of  the  pulse 

void 

void  space  of  region  filled  with  electrolyte 

total 

total  contributions  to  cell  sandwich  ASI 

+ 

lithium  ions  in  the  electrolyte 

s  lithium  in  the  active  material  particles 

T  maximum  lithium  in  the  active  material  particles 


a  1  cm2  section  of  electrode  to  be  constant  for  a  specific  chemistry 
and  state  of  charge  [4-8].  The  ASI  for  the  required  duration  of  the 
pulse  could  then  be  accounted  for  by  selecting  the  appropriate  data 
point  from  the  experimental  measurement  or  using  an  equivalent 
circuit  model  to  fit  the  experimental  results.  The  selected  value 
was  then  added  to  a  calculated  current  collection  system  resistance 
based  on  a  chosen  cell  design.  While  this  approach  is  fast  and  rea¬ 
sonable  under  certain  conditions,  the  assumption  of  a  constant  ASI 
for  many  disparate  battery  designs  is  physically  unrealistic  and  may 
result  in  the  prediction  of  batteries  that  are  impossible  to  reproduce 
with  experiments.  The  goal  of  this  work  is  to  present  an  improved 
approach  to  treating  the  ASI  that  maintains  the  speed  and  simplic¬ 
ity  required  while  incorporating  the  necessary  limitations  resulting 
from  the  underlying  electrochemical  processes.  We  will  use  Li-ion 
chemistries  as  the  basis  for  discussion  throughout  the  remainder 
of  the  paper;  although,  the  approach  we  will  present  is  sufficiently 
general  to  be  applied  to  a  range  of  battery  chemistries. 

2.  Model  development 

2A.  Formulation  of  ASI  calculation 

A  schematic  of  a  section  of  a  Li-ion  cell  is  depicted  in  Fig.  2.  Ah 
three  regions  contain  a  void  or  porous  phase  that  is  filled  with  elec¬ 
trolyte.  The  separator  region  contains  a  polymer  phase  to  prevent 
electron  transport  between  the  positive  and  negative  electrodes 
while  allowing  electrolyte  transport  between  these  two  regions. 
Each  electrode  is  composed  of  solid  particles  of  active  material 
for  intercalation  and  inactive  materials  that  improves  electronic 
conduction  and/or  structural  integrity.  The  active  materials  are  typ¬ 
ically  covered  in  a  thin  him  formed  from  side  reactions  with  the 
electrolyte  referred  to  as  the  solid  electrolyte  interphase  (SEI).  The 
electrons  move  to/from  the  electrodes  to  their  respective  current 
collectors  and  in/out  of  the  cell  to  the  external  circuit. 


Fig.  1.  Controlled  current  and  observed  voltage  during  a  HPPC  test. 
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Fig.  2.  An  illustration  of  a  segment  of  a  full  Li-ion  cell.  The  volume  fractions  in  each 
region  sum  to  unity.  All  regions  contain  a  void  or  porous  phase  that  is  filled  with 
electrolyte. 


The  ASI  is  an  experimentally  measured  value  resulting  from  a 
collection  of  resistances  that  originate  from  transport  and  kinetic 
processes  within  the  current  collection  components,  electrolyte, 
active  material,  and  their  shared  interfaces.  The  contribution  to  the 
ASI  from  the  current  collection  system  should  be  calculated  sepa¬ 
rately  to  maintain  simplicity  and  ease  of  use  of  the  model.  A  large 
portion  of  the  ohmic  contribution  to  the  ASI  arises  from  ion  trans¬ 
port  through  the  electrolyte  filled  porosity  of  the  electrode  and 
separator,  and  any  SEI  that  may  be  present  on  the  active  materi¬ 
als.  An  approximation  of  these  resistances  is  found  in  Eq.  (2)  for  the 
cell  sandwich  depicted  in  Fig.  2.  This  equation  is  only  to  serve  as  an 
introduction  to  some  of  the  straightforward  sources  of  impedance. 
The  region  thickness,  L,  void  volume  fraction,  £void,  electrochemi- 
cally  active  area  per  unit  volume,  a,  and  thickness  of  the  SEI,  8sei , 
are  most  likely  different  for  the  positive  and  negative  electrodes. 
The  conductivity  of  the  bulk  electrolyte,  /c0,  is  significantly  greater 
than  the  SEI,  Kse j.  The  ASI  from  ion  transport  within  the  pores  of  the 
electrode  scales  with  one  half  of  the  electrode  thickness  assum¬ 
ing  a  constant  current  distribution  throughout  the  depth  of  the 
electrode.  Ohmic  losses  resulting  from  the  transport  of  electrons 
through  the  electrode  conducting  matrix  to  the  current  collector 
must  also  be  included  in  the  ASI.  The  treatment  of  those  will  be 
discussed  later  in  the  paper. 


ASIU+  < 


4m 


K™aL 


2^oid 


Jsei_ 
Ksei  al 


(2) 


The  charge  transfer  resistance  (CTR)  component  of  the  ASI  is 
perhaps  the  most  variant  with  changes  in  electrode  design.  We 
may  examine  the  CTR  by  analyzing  the  Butler-Volmer  equation 
often  used  to  predict  electrochemical  reaction  rates.  The  nonlin¬ 
ear  dependence  of  current  on  potential  can  be  simplified  to  a  linear 
approximation  for  this  discussion.  An  example  of  this  linear  formu¬ 
lation  of  electrochemical  reaction  rate  per  unit  of  interfacial  area, 
in,  is  expressed  in  Eq.  (3)  [9]. 


in  =  io 


/  cT  -  cs  \  /  cs  \“c  (aa  +  ac)F 

Ut-cs*;  Uv  rt  ^ 


(3) 


Here,  i0  is  the  exchange  current  density  related  to  the  interfacial 
area.  R  and  T  correspond  to  the  universal  gas  constant  and  abso¬ 
lute  temperature  respectively.  The  subscripts  on  the  concentration 
terms,  Cj ,  are  as  follows:  +  is  for  the  lithium  ions  in  the  electrolyte,  s 
is  for  lithium  in  the  solid  active  material  and  T  is  for  the  maximum 
lithium  content  in  the  solid.  The  surface  overpotential,  r /,  is  the 


difference  between  the  electrode  potential,  V,  and  the  equilibrium 
potential,  U.  The  transfer  coefficients,  aa  and  ac,  are  each  taken  to 
be  0.5  and  thus  sum  to  unity. 

One  must  consider  the  relationship  of  the  transfer  current  den¬ 
sity,  in,  to  the  porous  electrode  of  the  electrochemical  cell  to  analyze 
the  CTR’s  role  in  the  ASI.  We  consider  the  positive  electrode  for 
the  first  example.  Porous  electrode  theory  relates  the  divergence 
of  the  superficial  electrolyte  current  of  an  electrode,  is,  to  the  rate 
of  electrochemical  reaction  in  Eq.  (4)  [10].  The  reaction  is  treated 
as  homogenous  by  the  inclusion  of  a,  the  ratio  of  electrochemi¬ 
cal  surface  area  to  electrode  volume.  The  negative  sign  in  Eq.  (4) 
reflects  the  arbitrary  decision  to  denote  a  discharge  current  as 
being  positive  valued.  We  choose  the  simplifying  assumptions  of 
a  unidimensional,  isothermal  model  while  also  neglecting  gradi¬ 
ents  within  the  electrode  to  allow  for  lumped  behavior.  That  is,  we 
can  consider  the  reaction  to  occur  equally  throughout  the  depth 
of  the  electrode,  Eq.  (5).  This  assumption  is  reasonable  for  small 
to  moderate  electrode  thicknesses  with  slow  interfacial  kinetics 
or  for  a  significant  sloping  open  circuit  potential  versus  state  of 
charge  curve  that  is  common  for  the  layered  oxide  intercalation 
compounds  commonly  used  in  Li-ion  batteries. 


The  current  density,  /,  is  related  to  the  geometric  area  of  the 
electrode  and  chosen  to  be  a  positive  value  for  discharge  currents 
and  negative  value  for  charge  currents.  The  concentration  of  the 
reactants  is  treated  by  accounting  for  the  depletion  of  the  species  at 
higher  current  densities.  Linear  concentration  profiles  are  assumed 
for  Li+  in  the  separator  and  within  the  solid  particles.  The  limiting 
current,  /lim,  for  lithium  transport  is  treated  with  a  simplification 
of  the  concentration  profile  as  a  function  of  current,  Eq.  (6). 


This  treatment  will  adequately  provide  for  a  dramatic  increase  in 
the  CTR  at  the  limiting  current  but  not  the  more  gradual  increase 
due  to  moderate  concentration  polarization.  Substitution  of  Eq.  (6) 
into  (5)  results  in  Eq.  (7). 


/  —  —  aioLpos  < 


(7) 


The  limiting  current  for  ionic  transport  in  the  electrolyte  and 
for  solid  state  diffusion  in  the  active  material  may  be  determined 
by  separate  calculations  using  diffusion  coefficients  and  diffusion 
lengths.  These  values  may  also  be  verified  with  a  detailed  model 
based  on  concentrated  solution  theory  and/or  with  experimental 
measurements.  Our  implementation  for  limiting  currents  repre¬ 
sents  a  best  case  scenario  and  most  likely  overestimates  the  limiting 
currents  for  transport  within  the  electrolyte  and  solid  active  mate¬ 
rial.  The  maximum  flux  of  Li+  through  the  solvent  filled  pores  of 
the  separator  may  be  approximated  using  dilute  solution  theory 
and  the  Bruggeman  correction  for  tortuosity  [10].  Eq.  (8)  uses  Fara¬ 
day’s  law  to  convert  the  maximum  flux  into  a  current.  Here,  DLiPF6 
is  the  diffusion  coefficient  of  a  commonly  used  electrolyte,  z+  is  the 
charge  of  the  cation,  and  t+  is  the  cation  transference  number  with 
respect  to  the  solvent.  The  change  in  cation  concentration  is  taken 
to  be  twice  the  initial  concentration,  Ac+  =  2c+i0.  The  role  of  elec¬ 
trolyte  transport  within  the  pores  of  the  electrode  is  not  treated  in 
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this  paper. 

, ionic  _  Z+FDLiPF6£TOUi  A C+ 

Mim  _  i  f  7  W 

The  analogous  limiting  current  for  lithium  diffusion  within  an 
intercalation  particle  is  Eq.  (9).  The  change  in  intercalated  lithium 
concentration,  Acs  =  Ct(1  -y),  and  solid  state  diffusion  coefficient, 
Ds,  both  depend  upon  the  normalized  amount  of  lithium  within 
the  particle,  y.  If  the  particle  is  assumed  to  be  at  50%  state  of  charge 
(SOC),  then  y  =  0.5.  This  equation  is  most  likely  a  poor  represen¬ 
tation  of  the  actual  limiting  current  as  the  boundary  conditions 
for  the  insertion  particle  are  different  than  those  of  transport 
through  the  porous  separator.  The  concentration  profile  within  the 
particle  is  most  likely  nonlinear  and  thus  an  alternate  approach  to 
determining  the  limiting  current  is  advisable.  We  choose  to  fit  our 
limiting  current  to  available  published  experimental  data. 

Afm iff  =  aLposZ+FDs  (9) 

rp 

The  electrochemical  active  area  per  unit  volume  may  be  related 
to  the  particle  radius  with  the  approximation  a  =  3eact/tp.  The  ben¬ 
efit  of  reducing  particle  size  is  two-fold.  The  smaller  particle  size 
decreases  the  time  constant  for  diffusion,  (rp)2/Ds.  Additionally,  the 
electrochemical  active  area  is  inversely  proportional  to  the  parti¬ 
cle  radius.  Thus,  a  smaller  particle  has  significantly  larger  surface 
area  for  the  electrochemical  reaction  to  occur.  The  practical  con¬ 
sequence  of  using  smaller  particles  is  that  the  electrode  becomes 
less  energy  dense.  The  electrode  void  fraction  is  always  higher  for 
an  electrode  composed  of  small  particles  rather  than  large  particles 
using  current  manufacturing  techniques.  A  balance  must  be  made 
between  the  rate  capability  of  an  electrode  and  the  desired  energy 
density. 

The  ASI  for  the  positive  electrode  CTR  is  related  to  the  kinetic 
expression  by  accounting  for  the  change  in  the  equilibrium  poten¬ 
tial  as  charge  is  passed,  Eq.  (10).  The  variables  on  the  right  hand  side 
of  the  following  equations  pertain  to  the  positive  electrode  unless 
otherwise  stated.  The  superscript  identifier  is  left  off  to  improve 
the  ease  of  reading.  The  concentration  of  lithium  at  the  surface  of 
the  particle  determines  the  equilibrium  potential.  One  must  know 
the  dependence  of  the  potential  on  the  SOC  and  how  the  concen¬ 
tration  of  lithium  in  the  particle  will  change  during  the  pulse  to 
calculate  the  thermodynamic  contribution  to  the  ASI.  The  equilib¬ 
rium  potential  function  for  a  positive  electrode  material  is  generally 
known  and  often  reported  in  literature.  We  will  use  the  slope  of  the 
known  potential  function  and  assume  it  is  constant  for  small  per¬ 
turbations  from  a  specified  SOC.  For  example,  a  10C  pulse  lasting 
1 0  s  is  assumed  to  be  a  small  perturbation  as  it  removes  less  than  3% 
of  the  total  lithium  of  an  electrode.  We  note  that  it  is  common  for 
potential  functions  to  rapidly  change  near  the  end  of  discharge.  This 
steep  change  will  manifest  as  a  large  contribution  to  the  ASI.  The 
example  we  explore  in  later  sections  is  based  on  materials  at  50% 
SOC  and  thus  the  slope  is  not  a  strong  function  of  concentration. 

...pcs  v^-vr  UtoOS-^°S  n  A  V  r,  dUdy 

ctr  1  I  I  I  1  dy  1 

(10) 

The  surface  concentration  of  the  particle  is  difficult  to  imple¬ 
ment  with  the  simplified  approach  used  in  this  work.  However, 
the  change  in  the  equilibrium  potential  with  the  average  concen¬ 
tration,  that  which  exists  after  a  long  relaxation  period,  is  easily 
determined.  To  present  the  flux  of  lithium  in  the  solid  in  a  more 
intuitive  way,  the  current  density,  /,  is  normalized  to  a  C-rate,  rc, 
using  Eq.  (11).  Here,  the  specific  capacity,  Q,  the  active  material 
density,  p,  active  material  volume  fraction,  sac t,  and  the  electrode 
thickness,  L,  are  used.  The  change  in  the  average  concentration  of 


lithium  that  arises  from  the  passage  of  current  is  shown  Eq.  (12) 
where  tpulse  is  the  length  of  the  current  pulse.  The  final  expression 
for  the  positive  electrode  ASIctr  is  represented  by  Eq.  (13).  Here, 
we  account  for  the  concentration  polarization  in  the  second  term, 
rill  while  accounting  for  the  thermodynamic  factor  with  the  third 
term,  dU/I. 


I  =  rcQpsactLp 

os 


dy  «  Ay 


rcf pulse 

3600 


ASIcpt°rs  = 


RT 

aioLpOSF 


d  U  fpulse 
dy  3600Qp£actf'pos 


1 


(11) 

(12) 


(13) 


In  practice,  it  is  difficult  to  separate  the  CTR  from  the  resistance 
of  the  transporting  ions  through  the  SEE  Both  scale  with  (aL)-1 
and  would  require  studies  such  as  electrochemical  impedance 
spectroscopy  (EIS)  to  separate  out  the  resistances  based  upon  the 
relative  time  scales.  Unfortunately,  EIS  does  not  unequivocally  pro¬ 
vide  an  answer  to  this  question  since  the  analysis  of  the  spectrum 
relies  upon  assumptions  of  which  processes  are  occurring.  We 
choose  to  sum  the  two  resistances  together  in  an  effort  to  accu¬ 
rately  capture  the  observed  behavior.  The  interfacial  ASI  for  the 
positive  electrode  may  then  be  considered  in  Eq.  (14). 


dU  t  pulse 

dy  3600Qp£act 


(14) 


The  remaining  contributions  to  the  ASI  must  be  accounted  for  to 
compare  Eq.  ( 14)  with  experimental  ASI  measurements.  We  choose 
to  lump  the  remaining  resistances  into  a  parameter,  ASICOnst,  which 
will  vary  only  with  differing  SOC.  Therefore,  all  changes  to  the  ASI 
with  electrode  thickness,  current  density,  or  C-rate  will  result  from 
the  interfacial  impedance  in  Eq.  (15). 

ASItotal  =  ASlf™  +  ASI^f  +  ASIconst  (15) 


2.2.  Determination  of  the  parameters 

The  parameters  used  in  Eq.  (15)  are  a  combination  of  measured, 
fitted,  and  calculated  values.  The  determination  of  interfacial 
parameters  such  as  resistivity  of  the  SEI  or  exchange  currents  is 
challenging.  Not  surprisingly,  the  reported  values  for  exchange 
currents  and  SEI  parameters  vary  widely  throughout  the  literature. 
Consequently,  we  choose  to  assume  one  of  the  sources  of  the 
resistance  is  dominate.  The  impedance  of  the  positive  electrode 
SEI  is  usually  small,  especially  at  the  beginning  of  life.  Therefore, 
we  neglect  the  SEI  contribution  to  the  positive  electrode  interfacial 
ASI  for  the  following  calculations  in  this  paper.  The  equilibrium 
potential  function  for  the  positive  electrode  is  taken  from  the  work 
of  Dees  et  al.  [11]  with  50%  SOC  is  defined  at  a  open  circuit  cell 
voltage  of  3.68  V.  The  slope  of  the  potential  function  is  approxi¬ 
mately  -1  at  this  SOC.  The  positive  electrode  exchange  current 
will  be  fit  to  experimental  ASI  measurements  and  thus  include  the 
contribution  from  the  SEI.  We  adopt  a  value  of  50  £2  cm2  for  the 
negative  electrode  interfacial  ASI  based  upon  the  single  particle 
measurements  of  Dokko  et  al.  [12-14].  This  impedance  is  based 
on  particle  surface  area  calculated  with  the  radius.  Therefore,  we 
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Electrode  thickness  /  pm 


Fig.  3.  Experimentally  measured  ASI  values  (open  circles)  from  [15]  compared  to  fit 
(solid  line)  from  Eq.  (15). 


use  the  particle  radius  for  our  chosen  graphite  to  estimate  the 
electrochemical  active  area.  We  ignore  limiting  currents  in  the 
negative  electrode  based  upon  the  focus  on  discharge  behavior 
and  high  performance  communicated  by  Dokko  et  al.  [13].  The 
diffusion  coefficients  for  the  electrolyte  are  taken  from  Dees  et  al. 
[11  ]  with  DLiPF6  =  1.1  x  10-6  cm2  s-1  for  the  electrolyte  used  in  the 
experimental  measurements.  We  use  the  experimental  data  for 
current  pulses  from  Abraham  et  al.  to  define  the  limiting  C-rate  as 
27b-1  [16]. 

Fig.  3  compares  the  fit  of  Eq.  ( 1 5 )  to  experimentally  measured  ASI 
values  for  2032-type  coin  cell  configurations  with  varying  positive 
electrode  thicknesses  from  the  work  of  Lu  et  al.  [15].  The  ASI  shown 
is  for  10  s  1.8C  discharge  pulses  from  a  coin  cell  at  50%  SOC.  The 
parameters  used  for  the  fit  are  shown  in  Table  1.  All  cells  are  based 
on  a  positive  electrode  active  material  of  LiNio.80Coo.10Alo.05O2 
(NCA)  layered  oxide  (Fuji  CA1505)  and  a  negative  electrode  based 
on  G8  graphite  from  Conoco  Philips.  The  manufacturer  reports  an 
average  particle  size  of  8  pm  and  a  specific  surface  area  of  2  m2  g-1 
using  the  BET  method.  The  thickness  of  the  graphite  is  varied  to 
maintain  an  excess  negative  electrode  capacity  of  1.1 -1.5  com¬ 
pared  to  the  positive  electrode.  The  separator  is  Celgard  2325  and 
the  electrolyte  is  1.2  M  LiPFg  in  EC: EMC  (3:7  by  weight). 

Good  agreement  is  achieved  between  the  model  and  the  full 
cell  data  with  only  two  fitting  parameters,  i0  and  ASICOnst.  The  ASI 
is  nearly  constant  for  electrode  thicknesses  greater  than  40  pm.  In 
contrast,  the  positive  electrode  dominates  the  ASI  at  small  electrode 
thicknesses  for  the  full  cell  measurements,  Fig.  4.  The  change  in 
positive  electrode  equilibrium  potential  is  responsible  for  approx- 


Table  1 

Parameters  used  in  Eq.  (15)  to  fit  experimental  data. 


Property 

Symbol 

Graphite 

NCA 

Units 

Specific  reversible  capacity 

a 

310 

150 

mAhg-1 

Exchange  current 

io 

N/A 

0.63 

mAcm2 

Particle  radius 

rP 

4 

5 

|jim 

Electrochemical  active  area 

a 

4667 

2893 

cm-1 

Active  material  density 

P 

2.25 

4.74 

gem-3 

Active  volume  fraction 

£act 

0.62 

0.48 

Void  volume  fraction 

£void 

0.32 

0.32 

Fig.  4.  The  total  ASI  is  composed  of  contributions  from  the  positive  and  negative 
interfacial  impedance  and  the  lumped  constant  parameter.  The  positive  electrode 
interfacial  impedance  dominates  the  total  ASI  at  low  electrode  thicknesses. 

imately  one  third  of  the  positive  electrode  interfacial  impedance. 
The  contribution  of  the  negative  electrode  interfacial  impedance  is 
never  more  than  25%  of  the  positive.  The  lumped  parameter  most 
likely  captures  some  of  the  complicated  interfacial  behavior  from 
both  electrodes. 

Abraham  et  al.  used  a  reference  electrode  cell  to  study  the  ASI  as 
a  function  of  pulse  magnitude  for  a  NCA/Mag-10  graphite  cell  [16]. 
Fig.  5  compares  the  measured  and  calculated  ASIs  for  1 8  s  discharge 
pulses  from  a  full-cell  balanced  to  a  35  pm  thick  positive  electrode 
at  50%  SOC.  The  limiting  C-rate  was  fit  to  the  experimental  data 


Fig.  5.  Comparison  of  the  ASI  for  18  s  discharge  pulses  from  [16]  (open  circles)  to 
the  model  (solid  line).  The  limiting  current  was  fit  to  the  values  seen  here. 
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shown  and  is  only  significant  close  to  the  limiting  value  of  27  h-1 . 
The  lumped  parameter,  ASICOnst,  was  increased  to  match  the  larger 
ohmic  losses  in  the  reference  electrode  cell  of  the  Abraham  et  al.  The 
larger  ohmic  losses  are  the  result  of  electrical  contact  resistances 
between  the  electrode  foil  and  current  collector  and  the  double 
separator  surrounding  the  reference  electrode.  The  good  fit  is  verifi¬ 
cation  that  the  approach  taken  in  this  model  is  sufficient  to  provide 
reasonably  accurate  results  to  feed  into  a  battery  design  model. 

2.3.  Limitations  of  the  model 

A  more  rigorous  approach  would  treat  the  competing  effects 
that  occur  within  the  porous  electrode.  The  current  distribution  in 
the  porous  electrode  is  rarely  ever  totally  uniform.  In  most  cur¬ 
rent  Li-ion  cells,  the  solid  matrix  conductivity  is  much  higher  than 
the  electrolyte  conductivity.  The  ease  of  electron  transport  in  the 
solid  phase  (i.e.  near  constant  potential)  has  the  effect  of  shift¬ 
ing  the  current  density  towards  the  separator-electrode  interface 
to  minimize  the  ohmic  losses  in  the  electrolyte.  The  set  of  equa¬ 
tions  used  to  accurately  predict  the  current  potential  behavior  has 
been  widely  accepted;  however,  many  parameters  require  deter¬ 
mination  and  computation  is  not  trivial  [9,10,17,18].  Euler  and 
Nonnenmacher  and  then  Newman  and  Tobias  treated  this  complex 
problem  with  simplifications  [10,19,20].  An  analytical  solution  has 
been  presented  that  allows  the  current  distribution  to  vary  along 
the  electrode  thickness  as  dictated  by  the  current  density,  reaction 
kinetics,  and  ohmic  losses.  To  achieve  this  end,  the  concentration 
of  the  reactant  is  assumed  constant  and  film  resistances  are  not 
treated.  We  have  compared  this  analytical  solution  with  our  more 
empirical  approach.  Our  proposed  ASI  calculation  achieves  a  better 
fit  with  the  experimental  data  than  the  analytical  solution  of  Euler 
and  Nonnenmacher  especially  in  the  transition  from  a  constant  ASI 
to  the  steep  increase  at  small  electrode  thicknesses.  This  agreement 
may  be  fortuitous,  but  the  failure  of  the  early  analytical  porous  elec¬ 
trode  solution  most  likely  results  from  the  complex  mass  transport 
occurring  within  the  pores  of  the  electrode,  intercalation  particles, 
and  the  film  that  serves  as  the  interphase. 

The  limitations  of  the  proposed  approach  are  a  result  of  the 
simplifications  used  to  enable  quick  calculations.  The  equation  is 
derived  for  the  linear  regime  of  the  Butler-Volmer  equation.  This 
typically  applies  for  low  to  moderate  C-rates  in  the  Li-ion  bat¬ 
tery  system.  The  exponential  dependence  of  the  full  Butler-Volmer 
equation  is  likely  to  be  important  for  higher  C-rates.  Furthermore, 
this  approach  assumes  an  isothermal  cell  design.  This  assump¬ 
tion  may  be  reasonable  if  the  designed  deviation  from  open  circuit 
voltage  is  a  low  value.  For  example,  we  typically  design  batter¬ 
ies  that  achieve  maximum  power  at  only  80%  of  the  open  circuit 
voltage.  Inclusion  of  the  full  Butler-Volmer  equation  and  the  heat 
generation  will  lower  measured  ASI  values.  However,  we  also  do 
not  account  for  uneven  current  distributions  in  the  porous  elec¬ 
trode  structure  that  will  increase  measured  ASIs.  In  addition,  the 
ionic  limiting  current  is  approximated  for  a  steady  state  condition. 
Higher  currents  are  most  likely  possible  for  short  durations  as  the 
concentration  of  the  reactants  will  be  higher  in  the  initial  moments 
the  current  is  passed.  The  level  of  accuracy  achieved  with  our  stated 
assumptions  is  reasonable  to  achieve  the  goal  of  the  model.  The 
intention  of  this  work  is  to  provide  a  realistic  result  that  can  guide 
the  battery  designer’s  intuition  or  feed  into  other  models. 

3.  Results  and  discussion 

3.1.  The  dependence  of  the  ASI  on  electrode  thickness  and  C-rate 


Fig.  6.  The  effect  of  limiting  currents  on  the  ASI  may  be  examined  by  changing  the 
C-rate  as  a  function  of  electrode  thickness. 

plays  the  change  in  calculated  ASI  as  the  C-rate  is  varied  from  5C 
to  20C.  The  two  different  limiting  currents  affect  the  ASI  in  dis¬ 
tinct  ways.  The  ASI  quickly  increases  to  an  infinite  value  as  the 
current  density  reaches  the  ionic  limiting  current  (i.e.  at  large  elec¬ 
trode  thicknesses).  The  maximum  achievable  electrode  thickness 
decreases  with  increasing  C-rate.  This  is  because  the  current  density 
passing  through  the  separator  increases  with  increasing  electrode 
thickness  (capacity),  at  constant  C-rate,  until  it  approaches  the 
limiting  current.  In  contrast,  the  ASI  at  all  electrode  thicknesses 
increases  as  the  C-rate  approaches  the  solid  state  diffusion  limited 
current.  In  this  analysis,  one  can  see  the  interplay  between  C-rate 
and  current  density.  In  actual  Li-ion  cells,  the  increase  in  impedance 
will  begin  at  much  lower  current  densities  and  C-rates  as  the  con¬ 
centration  gradients  build  in  the  cell.  The  goal  of  our  calculation 
is  to  set  a  physical  limitation  to  the  cells  one  might  design,  not  to 
precisely  capture  the  current-voltage  behavior. 

An  interesting  feature  of  the  ASI  at  constant  C-rate  is  the  increase 
in  impedance  as  the  electrode  thickness  decreases.  Dees  et  al.  pre¬ 
viously  simulated  the  same  NCA-graphite  system  with  a  robust 
model  based  on  porous  electrode  and  concentrated  solution  theo¬ 
ries  [11].  Their  results  also  predict  the  steep  increase  in  ASI  with 
decreasing  electrode  thickness.  The  work  presented  here  explains 
the  consequences  of  this  increased  ASI  in  the  context  of  the  end 
battery  design.  This  behavior  is  a  direct  result  of  the  inverse  rela¬ 
tionship  of  the  ASI  to  electrode  thickness  contained  in  Eq.  (14). 
At  constant  C-rate,  the  current  density  decreases  as  the  electrode 
thickness  decreases  and  thus  the  rise  in  ASI  is  cancelled  out  by  the 
diminishing  current.  This  is  best  demonstrated  by  examining  the 
voltage  loss  incurred  from  the  positive  electrode  CTR.  Multiplica¬ 
tion  of  Eq.  (13)  by  (12)  results  in  Eq.  (16)  with  f[I)  representing 
the  limiting  current  expression.  Inspection  of  Eq.  (16)  reveals  the 
voltage  loss  from  the  interfacial  impedance  does  not  depend  on 
the  electrode  thickness  if  the  C-rate  is  held  constant.  Clearly,  volt¬ 
age  losses  will  increase  if  the  current  density  is  held  constant  as 
electrode  thickness  is  reduced. 


The  effect  of  electrode  thickness  and  C-rate  on  the  ASI  may  be 
evaluated  using  Eq.  (15)  for  the  NCA-graphite  full  cell.  Fig.  6  dis- 


AV  =  rcQpsact 
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d  U  Fctpulse 
dy  3600 


(16) 
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Fig.  7.  The  effect  of  the  ASI  on  cell  voltage  at  the  end  of  a  10  s  IOC  discharge  pulse 
for  a  NCA-graphite  full  cell. 


Fig.  7  displays  the  cell  voltage  at  the  end  of  a  10C  pulse  that  is 
10  s  in  duration  as  a  function  of  positive  electrode  thickness.  This 
is  the  Vti  value  used  in  Eq.  (1).  The  cell  voltage  was  determined  by 
using  the  open  circuit  voltage,  Vt0  =  Ut0  =  3.68  V,  for  a  NCA-graphite 
full  cell  at  50%  SOC.  The  ASI  is  also  shown  for  the  same  conditions. 
The  steep  rise  in  the  ASI  at  low  electrode  thickness  does  not  result 
in  increased  polarization  of  the  cell  unlike  at  large  thicknesses.  On 
the  other  hand,  the  limiting  current  provides  for  a  severe  drop  in 
cell  voltage  at  the  maximum  electrode  thickness  for  this  C-rate. 

3.2.  Implementation  in  battery  design  model 


mined  by  the  specific  capacity  of  the  negative  electrode  and  the 
designed  excess  capacity  to  prevent  lithium  plating  during  charg¬ 
ing.  We  have  chosen  a  ratio  of  1 .25  negative  to  positive  capacity  for 
this  study. 

LP°S  =  QpSactApos  (18) 

We  have  calculated  a  number  of  high  P/E  ratio  batteries  to  high¬ 
light  the  role  of  changes  in  the  interfacial  ASI  on  battery  design. 
All  batteries  have  a  specified  power  of  25  kW  reached  at  80%  of  the 
open  circuit  voltage  using  one  insulated  module  containing  20  cells 
connected  in  series.  We  define  the  power  of  a  battery  based  upon 
the  pulse  power  characteristics  at  a  specific  SOC.  Batteries  used  in  a 
hybrid  electric  vehicle  require  a  high  power  to  energy  (P/E)  design 
and  typically  operate  at  50%  SOC.  The  usable  energy  of  these  cells 
is  assumed  to  be  25%  of  the  total  energy  measured  at  a  C/3  dis¬ 
charge  rate.  The  capacity  of  each  cell  changes  to  meet  the  stated 
energy  requirement.  Our  battery  design  model  explicitly  accounts 
for  the  physical  dimensions  and  weights  of  all  components  such 
as  cell  materials,  current  collection,  module  casing,  and  thermal 
insulation,  for  example.  The  results  for  a  NCA-graphite  system  are 
displayed  in  Fig.  8. 

The  FreedomCAR  energy  storage  goals  target  P/E  ratios  of  1 6-20 
for  power-assist  hybrids  batteries  [1].  For  the  battery  chemistry 
studied  here,  a  local  minimum  in  battery  weight  and  volume  is 
found  near  a  P/E  ratio  of  17h-1  (1.8  kWh).  The  maximum  power 
density  (i.e.  minimum  weight  and  volume)  is  located  at  the  knee 
in  the  ASI  curve.  The  battery  design  model  reduces  the  electrode 
thickness  as  the  required  energy  is  lowered.  The  decreasing  elec¬ 
trode  thickness  causes  the  ASI  to  increase  as  the  CTR  of  the  positive 
electrode  becomes  significant  as  discussed  earlier.  Fig.  9  illustrates 
the  dramatic  increase  in  cell  area  as  the  energy  requirement  and 
thus  electrode  thickness  are  reduced.  This  relation  is  a  direct  result 
of  the  interdependence  of  the  ASI,  L,  and  Acell  as  discussed  earlier. 

Fig.  10  clearly  shows  the  C-rate  is  the  only  significant  source  of 
polarization  from  the  two  limiting  current  mechanisms  considered. 
The  C-rate  increases  as  the  P/E  ratio  increases  at  lower  designed 
battery  energies.  The  current  density  remains  less  than  one  half  of 


The  role  of  changes  in  the  ASI  on  battery  design  may  be  eval¬ 
uated  using  the  Argonne  Battery  Design  model  which  has  been 
communicated  previously  [4-8  ].  The  fitted  ASI  equation  from  above 
is  implemented  into  an  iterative  scheme  that  calculates  current 
collection  resistance,  material  quantities,  and  battery  size  for  a 
specified  power  and  energy.  The  present  cell  design  is  based  on 
a  prismatic  pouch  cell  configuration  with  a  thickness  of  10  mm 
and  a  length  to  width  aspect  ratio  of  1.3.  The  positive  and  nega¬ 
tive  electrode  tabs  are  drawn  out  of  opposite  ends  and  are  nearly 
as  wide  as  the  cell  to  minimize  ohmic  losses.  The  current  collection 
components  depend  on  the  size  of  the  cells  but  typically  add  only 
~0.5  £2  cm2  to  increase  the  ASItotal  to  the  ASIbatt  for  the  batteries 
designed  below  (calculation  not  shown  for  brevity).  The  area  of  the 
positive  electrode,  Eq.  (17),  is  found  from  ASIbatt,  battery  power,  P, 
voltages,  V ,  and  number  of  cells,  Ncell.  The  difference  between  the 
open  circuit  voltage,  V0cv,  and  the  cell  voltage,  Vcen,  is  the  polariza¬ 
tion  of  the  cell.  The  negative  electrode  is  taken  to  be  1  mm  longer 
than  the  positive  electrode  in  both  length  and  width  to  alleviate 
concerns  of  lithium  plating  during  charge  pulses.  The  separator 
area  is  slightly  larger  than  the  negative  electrode  to  prevent  the 
electrical  shorting. 


_  _ ASIbatt;P _ 

P0S  “  Ncell\4ell(Vocv  -  Vcell) 


(17) 


The  positive  electrode  thickness,  Eq.  (18),  is  determined  from 
the  capacity  of  the  cell,  C,  specific  capacity  of  the  electrode,  and  the 
positive  electrode  area.  The  negative  electrode  thickness  is  deter¬ 


Energy  /  kWh 


Fig.  8.  Effects  of  battery  pack  energy  on  battery  volume  and  weight  for  20-cell, 
NCA-graphite  packs  delivering  25  kW  power  surges  for  10  s  at  80%  of  OCV. 
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Fig.  9.  Effects  of  battery  pack  energy  on  positive  electrode  area  and  thickness. 


the  limiting  value  and  is  dramatically  reduced  as  the  area  of  the  cell 
is  increased  (e.g.  ASI  increases).  In  fact,  the  total  current  is  the  same 
for  all  battery  designs  here  as  a  result  of  specifying  that  maximum 
power  will  be  achieved  at  a  set  fraction  of  the  open  circuit  voltage. 
Therefore,  doubling  the  number  of  cells  in  a  battery  from  20  to 
40  reduces  the  capacity  of  the  cell  by  half  through  a  reduction  in 
cell  area.  The  calculated  electrode  thickness  and  current  densities 
remain  the  same  as  they  are  governed  by  the  fraction  of  the  open 
circuit  voltage  at  which  the  maximum  power  is  achieved. 

The  effect  of  the  positive  electrode  CTR  on  battery  design  may 
be  examined  by  changing  various  parameters  such  as  the  solid 
state  diffusion  coefficient  and  the  particle  size.  Fig.  1 1  shows  the 


Fig.  10.  Current  density  (thin  grey  line)  and  C-rate  (thin  black  line)  as  a  function  of 
battery  energy  for  the  25  kW  design.  The  limiting  current  density  (thick  grey  line) 
and  limiting  C-rate  (thick  black  line)  are  also  shown  for  comparison. 


Energy  /  kWh 


Fig.  11.  The  effect  from  increasing  the  solid  state  diffusion  coefficient  (10DS)  and 
separately  reducing  particle  size  (10_1rp)  is  shown  on  the  battery  design.  Reducing 
particle  size  increases  the  limiting  C-rate  and  reduces  the  charge  transfer  resistance. 


change  in  battery  weight  and  volume  with  10DS  and  10-1  rp.  These 
choices  reflect  possible  material  changes  that  could  be  used  to  meet 
battery  power  and  energy  density  requirements.  The  reduction  in 
particle  size  has  the  strongest  effect  on  battery  design  by  lower¬ 
ing  the  impedance  and  thus  weight  for  all  P/E  designs  considered. 
The  change  in  void  volume  or  porosity  of  the  electrode  with  par¬ 
ticle  size  has  not  been  included.  This  could  be  easily  implemented 
with  an  empirical  correlation,  but  the  change  in  void  fraction  from 
33%  to  50%  has  only  a  minor  effect  on  battery  size  compared  to  the 
decrease  in  ASI.  Reducing  particle  size  may  result  in  other  detri¬ 
mental  effects.  The  increase  in  specific  surface  area  of  the  particle 
will  increase  side  reactions  with  the  electrolyte.  This  will  amplify 
the  rate  of  power  and  capacity  fade  while  possibly  exacerbating 
safety  issues.  Therefore,  blending  of  electrode  materials  has  been 
implemented  to  include  materials  that  have  high  power  density 
with  those  that  have  high  energy  density  [21  ]. 

The  operating  conditions  of  a  battery  system  should  be  chosen 
with  sufficient  room  for  changes  to  the  limiting  currents.  As  electro¬ 
chemical  systems  age,  the  porosity  of  the  separator  may  decrease 
and  some  fraction  of  the  intercalation  materials  may  become  inac¬ 
tive.  These  two  symptoms  of  aging  will  lower  the  limiting  current 
of  the  system.  The  performance  of  the  battery  will  be  severely  hin¬ 
dered  if  the  designed  operating  current  reaches  a  value  close  to  the 
limiting  currents.  The  design  of  an  electrochemical  system  should 
always  consider  the  lifetime  capability  of  the  system  rather  than 
only  initial  performance. 


4.  Conclusions 

A  simplified  expression  for  the  ASI  has  been  presented  and 
implemented  into  a  battery  design  model.  The  accuracy  of  pre¬ 
dictions  made  by  the  model  has  been  greatly  improved  by 
implementing  physical  limitations  in  the  interfacial  impedance 
expression.  This  iterative  approach  is  very  efficient  and  allows 
for  rapid  study  of  a  number  of  different  scenarios.  Furthermore, 
a  battery  design  model  with  this  simplicity  may  easily  feed  infor¬ 
mation  into  other  models  such  as  economic  plant  design  models 
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or  other  simulations  that  use  many  repetitive  calculations  [22,23]. 
The  combination  of  accuracy  and  ease  of  calculation  help  to  connect 
laboratory  studies  to  commercial  product  development. 
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